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ABSTRACT 

It has been proposed that the Hercules stream, a group of co-moving stars in 
the Solar neighborhood offset from the bulk of the velocity distribution, is the 
result of resonant interactions between stars in the outer disk and the Galactic 
bar. So far it has only been seen in the immediate Solar neighborhood, but the 
resonance model makes a prediction over a large fraction of the Galactic disk. I 
predict the distribution of stellar velocities and the changing Hercules feature in 
this distribution as a function of location in the Galactic disk in a simple model 
for the Galaxy and the bar that produces the observed Hercules stream. The 
Hercules feature is expected to be strong enough to be unambiguously detected 
in the distribution of line-of-sight velocities in selected directions. I identify quan- 
titatively the most promising lines of sight for detection in line-of-sight velocities 
using the KuUback-Leibler divergence between the predictions of the resonance 
model and an axisymmetric model; these directions are at 250° < / < 290°. The 
predictions presented here are only weakly affected by distance uncertainties, 
assumptions about the distribution function in the stellar disk, and the details 
of the Galactic potential including the effect of spiral structure. Caia and fu- 
ture spectroscopic surveys of the Galactic disk such as APOGEE and HERMES 
will be able to robustly test the origin of the Hercules stream and constrain the 
properties of the Galactic bar. 

Subject headings: Galaxy: bulge — Galaxy: disk — Galaxy: evolution — Galaxy: 
fundamental parameters — Galaxy: kinematics and dynamics — Galaxy: struc- 
ture 



1. Introduction 

The velocity distribution of nearby stars contains various large overdensities of co- 
moving stars that cannot be explained as the vestiges of clumpy star formation (e.g.. 
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Bovy & Hogg 2010). In particular, the Hercules moving group, a group of stars sig- 
nificantly offset from the bulk of the observed velocity distribution (e.g., Dehnen 1998; 
Bovy et al. 2009b, see Figure 1), displays a wide range of ages and metallicities (Blaauw 
1970; Raboud et al. 1998; Caloi et al. 1999; Bensby et al. 2007; Bovy & Hogg 2010). A 
compelling model has the Hercules stream originating through resonant interactions of 
stars in the outer disk with the bar in the central region of the Galaxy (Dehnen 2000; 
Fux 2001), but there is no confirmation of this picture beyond the locally observed stellar 
kinematics that it was proposed to explain. An investigation of the Hercules members' 
color-magnitude diagram found no evidence for any significant metallicity anomaly with 
respect to other local stars (Bovy & Hogg 2010), even though such a difference should be 
expected as the Hercules stars originate from a few kiloparsec toward the Galactic center 
where the average metallicity is higher than in the Solar neighborhood. The absence of the 
predicted metallicity anomaly could potentially be explained, e.g., through radial mixing 
(Sellwood & Binney 2002), but it is clear that chemo dynamical modeling of the Galactic 
disk is not currently up to the task of separating out these effects (e.g., Samland & Gerhard 
2003; Roskar et al. 2008). 

Here I propose that a robust and unambiguous test of the bar-resonance model for 
the Hercules stream is a search for the distinct pattern it predicts as we trace the velocity 
distribution of stars near the Solar radius around the Galaxy. Gaia is in a good position 
to have the final word on this, but I also show that the Hercules feature can be detected 
in the line-of-sight velocity distribution in selected directions on the sky. Clear signatures 
are predicted to be detectable in regions a few kiloparsec from the Sun, allowing for a clean 
verification of the bar-resonance model for the Hercules stream. 



2. Methodology 

To simulate the effect of the bar on stellar orbits near the Solar circle, we follow the 
approach of Dehnen (2000). For each fixed position in the Galaxy, we evaluate the velocity 
distribution function at values for the radial — toward the Galactic center — and tangential — 
in the direction of Galactic rotation — velocity by backward-integrating these velocities in 
the Galactic potential to obtain the intial orbit before bar formation. We posit that at this 
epoch the old stellar disk was in a steady-state that can be described by a simple distribution 
function that is a function of the energy and angular momentum of the orbit alone. It follows 
from the collisionless Boltzmann equation that the value of the velocity distribution at the 
present position and velocity is equal to that of the steady-state distribution function before 
bar formation evaluated at the initial position and velocity — or, equivalently, at the initial 
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orbit since the initial distribution function is time-independent. 

Since the Hercules stream is only apparent in the planar motions in the disk (e.g., 
Bovy et al. 2009b), we only consider a two-dimensional model for the Galaxy. We use a 
simple power-law rotation curve to model the axisymmetric Galactic potential 

v,{R) = vo{R/Rof , (1) 

where i?o is the distance from the Sun to the Galactic center. For most of the simulations 
below we specify this model further to a flat rotation curve (/3 = 0). The model for the bar 
is given by the following potential 



-{Rt/R)% for R>R^ 



w) = A.(t) com- m X I ,^;-;;3-'_2 ,Z'r<Z 

Here Qi^ is the pattern speed of the bar and i?^ is the bar radius, which is fixed to be 80 percent 
of the bar's corotation radius. The bar is grown smoothly during a time ti, always set to 
half of the total integration time in what follows, using the prescription 

A(t)^air^Vx|^^^-i^^ + i^ + ^' ^-'i-'^ (3) 

^ ^ 3 \Rj \ 1, t>h. ^ ^ 

A comparison of the local velocity distribution obtained by using this simple bar potential 
with that resulting from the more realistic one of Gardner & Flynn (2010) shows that this 
simple potential is sufficient to model the bar's influence on stellar orbits in the outer disk. 

For the distribution function of the stellar disk before bar formation, we use a Dehnen 
distribution function (Dehnen 1999a) given by 



fDehnen{E, L) OC ^ CXp 



S(i?e) 
4(^e) 



f^(i?e) [L - L,{E)] 



(4) 



4(^e) 

where i?e, Lc-> and Vt{Re) are the radius, angular momentum, and angular frequency, respec- 
tively, of the circular orbit with energy E. Using the procedure given in § 3.2 of Dehnen 
(1999a), we choose the and (Jr{R) functions such that they reproduce a disk with 

exponential surface density and velocity dispersion profiles 

= So exp {-R/Rs) , (7r[R) = ao exp {-R/R,) , (5) 

to an accuracy of a fraction of a percent at all radii. Our fiducial model for the distribution 
function has Rs= Ro/3 and i?cr= 3/?^, and ctq such that (Jr{Ro) = 0.2vo. 

As shown by Dehnen (1999b, 2000), the model described above with the Sun at Ro= 
0.9i?oLR — where i?oLR is the radius at which a circular orbit is in the outer Lindblad reso- 
nance with the bar — at a current angle of 25° with the bar, with a bar strength a = 0.01, 
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and integrating for four bar periods, reproduces the Hercules feature in the local velocity 
distribution while agreeing with other photometric and kinematical observations of the bar 
(e.g., Binney, Gerhard, & Spergel 1997; Bissantz & Gerhard 2002; Cole & Weinberg 2002; 
Shen et al. 2010). In what follows, unless indicated otherwise, we fix all of the parameters 
of the model at these fiducial values. All velocities are always considered to be with respect 
to their local standard of rest in the absence of a bar. 



3. The full planar velocity distribution in the disk 

Dehnen (2000) only simulated the velocity distribution at the position of the Sun for 
various combinations of the model parameters defined in the previous Section. Nevertheless, 
his results on how the velocity distribution at the Sun changes with the angle between 
the Sun-Galactic center line and the bar could also be interpreted as showing the velocity 
distribution on the Solar circle for diflFerent Galactocentric azimuths, fixing the bar angle. In 
this way, these simulations can make predictions for what we might find when we study the 
velocity distribution in regions beyond those probed by the Hipparcos mission (ESA 1997). 

In Figure 2, I present the predictions for the two-dimensional velocity distributions near 
the Solar circle for the full 360° range in Galactocentric azimuth for the fiducial model from 
§ 2. Only the range —90° < < 90° is shown here, since the barred Galaxy is symmetric 
with respect to rotations through 180°; the predictions for the other side of the Galaxy can 
be obtained from those shown by shifting the azimuths by 180°. The predicted velocity 
distributions on the Solar circle are not shown here, but they can be inferred from Figure 2 
in Dehnen (2000). 

It is clear from this Figure that we can trace the Hercules stream going around the 
Galaxy and when looking further in toward the Galactic center or further out. Furthermore, 
the effect of the bar is not confined to a narrow range in radii, but rather the effects of the 
bar and especially its outer Lindblad resonance are felt through a large range of radii (0.6 
< R/ Ro < 1.4). These effects are strong and robust predictions of the bar-origin model 
for the Hercules stream. This is especially the case in regions beyond the Solar radius. 
There the Hercules feature can be found significantly removed from the bulk of the velocity 
distribution, which might be severely shaped at all locations by the processes that create the 
pattern of low- velocity moving groups in the local velocity distribution (the main mode in 
Figure 1). As I discuss below, Gaia will be able to verify these predictions. 
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4. Line-of-sight velocity distributions 

Rather than wait for a comprehensive survey of the Galactic disk that measures the fuU 
phase-space position of a large number of stars by measuring proper motions, we can ask 
whether a targeted survey of line-of-sight velocities in selected regions of the Galaxy could 
detect a clear signature of the Hercules feature? The Hercules feature stands out in the 
velocity distribution most in its tangential velocity On the Solar circle near the Sun the 
line-of-sight velocity is almost entirely this Galactocentric tangential velocity so we would 
expect to see at least a hint of the Hercules stream in the line-of sight. 

In Figure 3 I show the predicted line-of-sight velocity distributions in various locations 
near the Solar circle all around the Galaxy obtained by integrating the full two-dimensional 
velocity distribution over the velocity component perpendicular to the line of sight. This 
Figure is no longer symmetric with respect to shifts of 180° in azimuth, since the line-of- 
sight direction breaks this symmetry. The gray curve is the line-of-sight distribution we 
would observe if the targeted stars were in a steady state in an axisymmetric potential 
with a Dehnen distribution function (equation [4]). Both line-of-sight velocity distributions 
are normalized to have equal mass since I assume that our knowledge of the Galaxy and 
the tracer population will not be good enough that we can test predicted differences in 
the absolute number of stars at a given location between the axisymmetric and the bar- 
resonance predictions. These differences are also largely a function of the particular form 
and parameters of the distribution function that we choose to represent the stellar disk before 
bar formation. 

The prospects to detect the Hercules feature are quite good on our side of the Galaxy. 
The Hercules feature is especially strong in the line-of-sight velocity distribution for Galactic 
longitudes 250° < / < 290°. This is because of the fortunate coincidence that for those lines- 
of-sight we see velocities in the direction along which the Hercules streams separates most 
from the bulk of the velocity distribution (see Figure 2) and that at these azimuths the 
Hercules stream is at the furthest distance from the rest of the velocity distribution. We see 
hardly any effect of the bar on the line-of-sight velocity distributions on the other side of the 
Galaxy (90° < (/) < 270°). On that side we mostly see the Galactocentric radial velocity in 
the line-of-sight, in which the Hercules feature does not separate clearly from the bulk of the 
stars. Even though upcoming surveys such as APOGEE have the ability to probe the other 
side of the Galaxy, for the purposes of detecting the Hercules stream this is uninteresting as 
it does not show up in the line-of-sight velocity distribution there. 

To further analyze the directions to target to best see the Hercules feature in the line-of- 
sight velocities, I have calculated the line-of-sight velocity distribution on a 100 by 20 grid in 
Galactocentric azimuth and radius. As a measure of the difference between the bar-resonance 
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model and the steady-state, axisymmetric model, I have calculated the KuUback-Leibler 
divergence between these predicted distributions (e.g., Mackay 2003) 



where, in this case, X IS the line-of-sight velocity, p{x) is the predicted distribution for the bar- 
origin model of the Hercules stream, and q{x) is the axisymmetric prediction. The KuUback- 
Leibler divergence is only one among many possible choices to quantify the difference between 
two probability distributions. It can be interpreted as the expected weight of evidence for 
the prediction by the bar-resonance model over the fully axisymmetric prediction that can 
be expected from each sampled star. As such, locations with the maximal Kullback-Leibler 
divergence are those in which the Hercules feature stands out most prominently in the line 
of sight. 

The result of this analysis is shown in Figure 4. I present results both using the full 
line-of-sight velocity distributions and only considering the part of the distribution at |'L'ios| > 
O.lbvQ. As is clear from Figure 1, the local velocity distribution contains many more features 
than the Hercules stream. These features are mostly contained to the main mode of the 
velocity distribution |'^| < 0.15'^o- Observations indicate that these features are also the 
result of dynamical interactions (e.g., Bovy & Hogg 2010; Sellwood 2010), albeit with spiral 
arms rather than the bar. Therefore, we can expect these or similar features to be present in 
the velocity distribution at other locations in the Galaxy as well, and thus, if the effect of the 
Hercules stream occurs close to the center of the bulk of the velocity distribution, it might 
be hard to disentangle from the other dynamical effects at play. However, these additional 
dynamical effects will be much less important for velocities > 0.15^o- Figure 4 shows that 
there are many locations, mostly at 250° < / < 290°, where the Hercules stream is clearly 
offset from the bulk of the velocity distribution and should be unambiguously detectable. 

Finally, in Figure 5 I show the influence of various nuisance and dynamical parameters 
on the predicted line-of-sight velocity distribution for the most promising region that is more 
than 3 kpc away. We see in this Figure that the effect of distance and realistic line-of-sight 
velocity uncertainties is negligible. The exact shape and parameters of the stellar distribution 
function also only minimally affect the predicted line-of-sight velocity distribution. Changing 
the scale length within its current uncertainty changes the prediction only in the amplitude 
of the effect: The Hercules stream is more prominent for shorter scale lengths since then 
there are relatively more stars in the inner Galaxy that are brought to the Solar circle by 
the interaction with the bar. Using a Shu distribution function (Shu 1969) instead of the 
Dehnen distribution function also gives very similar predictions. Changing the axisymmetric 
potential by adjusting the shape of the rotation curve shifts the peak as expected (see Dehnen 
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2000), but the Hercules feature nevertheless remains detectable. The properties of the bar 
itself, e.g., its pattern speed through i?oLR oi its strength a, change the predicted pattern, 
but they do not significantly alter the prospects for detection. 

5. Discussion 

The discussion in the previous sections shows that the bar-origin paradigm for the 
Hercules moving group makes strong predictions about how the Hercules feature should show 
up when we probe the velocity distribution in parts of the Galaxy far removed from the Solar 
neighborhood. In order to most strongly test the bar origin for the Hercules moving group 
it is best to extrapolate this picture far from the Solar neighborhood and look at regions 
further than a few kpc removed. Gam's parallaxes, proper motions, and mission- averaged 
line-of-sight velocities for dwarfs nearby and giants at distances greater than a few kpc will 
allow the full two-dimensional velocity distribution to be confronted with the predictions of 
Figure 2 (see below). However, the final Gaia catalog will only be available around 2020. 

The predicted line-of-sight velocity distributions in the previous section show that the 
bar-origin theory for the Hercules moving group can be tested with line-of-sight velocities 
alone. At distances greater than a few kpc this will mean observing giants to trace the 
velocity distribution. Large, planned spectroscopic surveys with line-of-sight velocity uncer- 
tainties less than 1 km s~^ can test the predictions made in this paper in the next few years 
(see below). These surveys will not need to observe many stars to test the predictions laid out 
here: Drawing samples from the predicted velocity distribution in Figure 5 shows that about 
500 stars suffice to distinguish the bar-origin model from the axisymmetric prediction (based 
on a KS test with P < 0.01). Since large distance uncertainties do not affect the prediction 
very much (see Figure 5), these stars can even be sampled throughout a volume of 500 pc or 
more. These spectroscopic surveys will also obtain detailed elemental abundances for their 
target stars, so they would also allow for a more detailed analysis of the abundances of stars 
in the Hercules stream (following, e.g., Bensby et al. 2007; Bovy & Hogg 2010), which could 
shed more light on radial mixing in the Galactic disk. 

The last two panels of Figure 5 show that the predicted line-of-sight velocity distribution 
depends strongly on the dynamical parameters of the bar. Changing the pattern speed of 
the bar through i?oLR shifts the predicted distribution and in particular the location of the 
Hercules feature. While this shift at a particular location is degenerate with the shift due to 
the uncertainty in the local circular velocity, results from various locations can be combined 
to break this degeneracy and the local circular velocity is also already quite well constrained 
(e.g., Bovy et al. 2009a). Changing the strength of the bar changes the relative heights of the 
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two main peaks in the predicted velocity distribution. However, the main peak of the velocity 
distribution might not be as simple as the one predicted here due to other dynamical effects. 
Integrated measures, such as the number of stars in the main peak and the number of stars 
in the Hercules feature, could mitigate this somewhat. While predicting actual constraints 
would involve detailed simulations of the expected data, it is clear that constraints on the 
dynamical properties of the bar can be derived from the experiment proposed in this paper. 

5.1. Gaia 

The astrometric Gaia mission will measure the parallaxes and proper motions of up to 
one billion stars, most of which will be disk stars (Ferryman et al. 2001). With parallax 
and proper motion accuracies down to 10 /xas and mission- averaged line-of-sight velocity 
uncertainties smaller than 10 km s~^ for many stars to 17th magnitude, Gaia will be able 
to probe the kinematics of the disk out to several kpc in all directions (Katz et al. 2004; 
Bailer- Jones 2008). Crucially, it will provide large samples of giants in many of the promising 
regions in Figure 2, with velocity measurements accurate enough to verify the predictions 
laid out in this paper. This will also allow the dynamical properties of the bar to be tightly 
constrained, something that will be hard for Gaia to do using direct observations of the 
bulge due to the large extinction toward the Galactic center (Robin et al. 2005). 

5.2. Radial velocity surveys 

Since the extinction in the optical bandpasses is rather large in the directions identified 
as promising in Figure 4 [Ay ^ 8 mag toward / ^ 270°; Marshall et al. 2006), spectroscopic 
surveys in the near-infrared such as AFOGEE are the best option for seeing the Hercules 
feature in the line-of-sight direction. The range of Galactic longitudes that AFOGEE can 
target along the disk is —5° < / < 250° and its expected line-of-sight uncertainties are about 
0.5 km s~^ (SDSS-III Froject Description-^). From Figure 4 this longitude limit is right on 
the edge of where we might hope to detect the Hercules feature. The Hercules bump in the 
direction / = 250° is about half as strong as in Figure 5, but remains as a clear signature of 
the bar. As before, about 500 stars a few kpc away are enough to detect the Hercules feature. 
Since AFOGEE is expected to start taking data in 2011 and to finish operations before the 
first Gaia data release, AFOGEE could potentially trace the Hercules stream around the 
Galaxy before Gaia does. 



Available at http://sdss3.org/collaboration/description.pdf . 
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Another planned spectroscopic survey is HERMES (Freeman et al. 2010). While oper- 
ating in the optical, it expects to survey hundreds of thousands of giants in the inner Galaxy, 
many as far as 6 kpc away, with line-of-sight velocity uncertainties less than 1 km s~^ (K. 
Freeman, private communication). Since HERMES operates from the Southern hemisphere 
it can observe the most promising directions in Figure 4 and, like APOGEE, detect the 
Hercules feature before Gaia. 

It is a pleasure to thank the anonymous referee for valuable comments and Mike Blanton 
and David W. Hogg for helpful discussions and assistance. Financial support for this project 
was provided by the National Aeronautics and Space Administration (grant NNX08AJ48G) 
and the National Science Foundation (grant AST-0908357). 
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Fig. 1. — Velocity distribution of nearby (< 100 pc) stars from Bovy et al. (2009b). vr and 
V(fy are the velocities toward the Galactic center and in the direction of Galactic rotation, 
respectively, corrected for the Solar motion. The Hercules moving group is visible as an 
over density at v^fy^ -0.2 '^o- 
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Fig. 2. — Two-dimensional velocity distributions along the Solar circle in the bar-resonance 
model for the Hercules stream. The x-axis in each subplot is the Galactocentric radial 
velocity and the ^-axis is the Galactocentric tangential velocity; the axis-ranges of each 
subplot are the same as in Figure 1. The velocity distributions on the other side of the 
Galaxy are obtained by shifting the azimuths by 180°. 
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Fig. 3. — Line-of-sight velocity distributions along the Solar circle. The range in line-of-sight 
velocities in each subpanel is as in Figure 5. The gray curve in each panel is the predicted 
distribution for an axisymmetric, steady-state Galaxy. 
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Fig. 4. — KuUback-Leibler divergence between the predicted line-of-sight velocity distribu- 
tions with and without the bar. Darker values have a larger discriminatory potential. The 
left panel uses the full line-of-sight velocity distributions, the right panel excludes the main 
mode of the velocity distribution that is likely to be affected by other dynamical and phase- 
space effects. The concentric dashed lines indicate constant-distance slices from the Sun, 
starting at 1 kpc and spaced 2 kpc apart. The white dashed lines indicate the boundaries 
of the four quadrants in Galactic longitude. The white dashed-dotted line shows the region 
of the sky that APOGEE can observe (-5° < / < 250°). 



- 15 - 




Fig. 5. — Variation of the predicted velocity distribution at i? = 1.075 i?o, = —21° {d= 0.4 
i?o, I = 270°) with distance uncertainties; hne-of-sight velocity uncertainties; parameters and 
shape of the distribution function; shape of the rotation curve; bar i?oLR; and bar strength. 



